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Abstract. The Orr-Sommerfeld equation is a spectral problem which is known to play an im- 
portant role in hydrodynamic stability. For an appropriate operator theoretical realization of the 
equation, we will determine the essential spectrum, and calculate an enclosure of the set of all eigen- 
values by elementary analytical means. 



1. Introduction 



The Orr-Sommerfeld equation is an non-selfadjoint general eigenvalue problem of 
the form 

(1.1) {-D 2 + a 2 ) 2 u + iaR[V ■ (-D 2 +a 2 )u + V" ■ u] = \{-D 2 + a 2 )u on I, 

subject to Dirichlet boundary conditions for u on the real interval I. Here, D := d/dx, i 
is the imaginary unit, and R > is the Reynolds number of an underlying fluid which 
moves in a stationary flow (perpendicular to /) with given real-valued flow profile 
V G C 2 (I). This flow is perturbed by a single-mode perturbation with wave number 
a > 0, and the physical question of stability or instability of the underlying flow in 
response to this perturbation arises. This question is closely related to the spectrum 
of (an appropriate operator-theoretical realization of) the Orr-Sommerfeld equation 
([L~l]). Essentially the flow is unstable (with respect to the wave number a) if (1.1) has 
an eigenvalue with negative real part. 



In this article, we will exclusively be concerned with the case I = [0, oo) correspond- 
ing to the half-plane flow along a wall (an overview o f res ults for this case can be found 
in ||). So the boundary conditions to be added to fll.l] ) read 

(1.2) u(0) = u'(0) = u(oo) = u'(oo) = 0. 



A flow profile V which is of particular interest in this case (but not under exclusive 
concern here) is the Blasius profile defined by V := /', where / is the solution of the 
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nonlinear boundary value problem (Blasius equation, a special case of the Falkner- 
Skan equation) 

(1.3) 2/"' + //" = on [0, oo), /(0) - /'(O) = 0, lim /'(*) = 1, 

x — >oo 

which can be shown to exist and to be unique, and moreover, to provide (for V = /') 

V — > 1, V — > 0, V" — > (exponentially) as x — > oo, 
y>0, V'>0 5 < on (0,oo). 



(1.4) 




Figure 1: Blasius profile together with it's first two derivatives 



We choose the following operator theoretical setting for problem (hi), (1.2), which 
we believe is the most natural and simple one: 

Let X and Y denote the complex Banach spaces H 2 (0, oo) P\Hq (0, oo) and L 2 (0, oo), 
respectively, and let D(A) := # 4 (0, oo)ni$(0, oo). The operators A : D(A) C X -> F 
and _B : X y are defined by Au and Bit denoting the left-hand and the right-hand 
side (without A) of equation ([hi]), respectively. It can be shown that A is closed (see 
III Theorem IX. 9. 5], plus the remark that the norm in X is stronger than the one in 
L 2 (0, oo)), and that B is one-to-one and onto, and bounded with bounded inverse. 
Our formulation of (|l . l|) , (1.2) now reads 



(1.5) 



Au 



\Bu. 



So we choose a direct operator theoretical realization of the Orr-Sommerfeld problem, 
rather than a formulation via an auxiliary operator constructed by Riesz's represen- 
tation lemma, as done e.g. in 

Our goal in this article is twofold: In Section 2, we will determine the essential 
spectrum of problem (|h^) exactly; the main tool is the well-known theorem on preser- 
vation of the essential spectrum under relatively compact perturbations. In Section 3, 
an enclosure for the set of all eigenvalues of problem (|l.5|) is calculated by elementary 
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analytical means. Here, exploiting sign restrictions on the profile V and its derivatives 
(such as flOp), we obtain results which extend enclosures known in the literature (see 
|2j, ||) , even if those have been obtained by more involved methods. In Q correspond- 
ing results have been obtained for the case of a compact interval /. For further results 
on the Orr-Sommerfeld equation on a compact interval see [|| and the literature cited 
there. 



2. The essential spectrum 



Since several definitions of essential spectra are around in the literature (see |l], 
Chapter IX] for an overview), we start with the definition which we are using. 

With X, Y denoting two complex Banach spaces, and A : D(A) C X — > Y a closed 
linear and B : X ^ Y a, bounded linear operator, we call the set 

<iess(A, B) := {A G C : A — XB is not a Fredholm operator of index 0} 

the essential spectrum of the problem Au = XBu or of the pencil (A, B). Here, a 
closed linear operator T : D(T) C X — > Y is called a Fredholm operator if its range is 
closed and has finite codimension d(T) in Y, and if its nullspace has finite dimension 
n(T); it has index if d(T) = n(T). Furthermore, 

g(A, B) :={A 6 C : A — XB is one-to-one and onto} 

denotes the resolvent set of the problem Au — XBu (or of the pencil(A, B)), and 
cr(A, B) := C \ g(A, B) its (total) spectrum. It is quite obvious that a ess (A,B) C 
a(A,B) and that a(A,B) \ aess(A,B) consists only of the eigenvalues of finite geo- 
metric multiplicity. 

We return to the specific choice for X, Y, A, B introduced in Section 1, i.e., describing 
the Orr-Sommerfeld problem. Our main result of the present section is 



Theore m 2 .1. Suppose that the flow profile V satisfies V — > 1, V" — > as x — > oo 
(co mpa re (|-?.^| )). With R and a denoting the parameters in (1.1), the essential spectrum 
of (|i. 4 ) i s given by 



(2.1) (Tess(A, B) = {/i + a 2 + laR ;/i€[0, oo)}. 

Proof. We introduce the auxiliary operators 

A : D(A) -> Y, A Q u := (-D 2 + a 2 ) 2 u + iaR(-D 2 + a 2 )u, 
K : D(A) Y, Ku := iaR[(V - 1) • (-D 2 + a 2 )u + V" ■ u] 

so that A = Ao + K. Since V — 1 and V" tend to zero and K is of lower order 
than Ao, it follows from [Q, Theorem IX. 8. 2] that K is relatively compact to Aq, i.e., 
for each || • ||x-bounded sequence (u n ) in D(A) such that (A u n ) is bounded in Y, 
(Ku n ) contains a convergent subsequence. Strictly speaking, the cited Theorem in ffl 
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provides the result with X and D(A) replaced by L 2 (0, oo) and H A (0, oo), respectively, 
but from this it is easy to obtain the same statement with X and D{A). 

Since B : X — > Y is one-to-one, onto, and bounded with bounded inverse, the above 
result implies that B~ X K : D(A) — » X is relatively compact to B~ 1 A$ : D(A) — > X. 
Now the invariance of the essential spectrum under relatively compact perturbations 
(see e.g. [Q, Theorem IX. 2.1]) provides 

(2.2) aessiB^Ao + B^K^dx) = aessiB^Ao^dx). 

The cited theorem in |l| is formulated for densely defined operators, but the proof 
clearly does not use this assumption. (This density is a general assumption in jl], 
Chapter IX] to be able to work with adjoint operators, which however do not occur 
in Theorem IX. 2.1 or in its proof.) Since B and B^ 1 are bounded, so that A — XB is 
a Fredholm operator of index if and only if B~ l A — Xidx is, we obtain furthermore 
aess(A, B) = (JcssiB^ 1 A, idx), so that (2.2) and the identity Ao + K = A provide 

(2.3) a eS s(AB) = aeas^AoMx)- 



It remains to be shown that tress {B 1 Ao, idx) equals the right-hand side of (2.1) which 
we call M from now on. This is done in two steps: We prove 

a) C \M C Q{B~ 1 A ,idx) (implying M D a(B- l Ao,idx) D 0-ess(fl -1 4>,*djc)). 

b) M C aessiB^Ao^dx). 

First we observe that, with E(x) := e~ ax , 

(2.4) B~ l A u = -u" + (a 2 + iaR)u + u"(0) ■ E for u G D(A), 

since the expression on the right-hand side is in X (for u G D(A)) and B applied to 
it equals Aqu. 

ad a) Let A G C \ M, so that A = fi + a 2 + iaR with /i G C \ [0,oo). Thus, 
(J, G g(C,id Y ) where C : D(C) C Y -> Y is given by D(C) := # 2 (0,oo) n f#(0,oo), 
Cu := -u". 

To prove that B~ x Aq — Xidx ■ D(A) — > X is onto, let r e X and define u := 
(C - Li)~ 1 r G D(C) = := (C - n)~ 1 E G D(C) = X. It is easy to calculate w in 

closed form and to show that w'(0) =^ 0. Therefore, u :— v — -jjftk w G X solves the 
equation 

, ^ „ v'(Q) 

(2.5) = r 

and moreover, u(0) = u'(0) = 0. Since u,E,r G H 2 (0, oo), ( |2.5| ) implies u" G 
# 2 (0,oo). Therefore u G ff 4 (0,oo). Altogether, it G D(A). Since u(0) = r(0) = 
and J5(0) = 1, (U) shows that u"(0) = u'(0)/u/(0), and (|J), Q then yield 
(S _1 i4o - A)u = r. 

To prove that B~ X A — Xidx is one-to-one, let u G -D(A) satisfy {B^ 1 A — X)u = 0, 
i.e., 

-u" - [iu + u"(0)E = 0. 
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Thus, u £ D(C) and (C-fi)u = -u"(0)E, i.e., u = -u"(0)w. Since u'(0) = 0, w'(0) + 
0, this implies u"(0) — and therefore u — 0. 

ad b) Let A £ M, so that A = /i + a 2 + iaR with li £ [0, co). Let 7 £ C and 
r := 7(1 — E)E £ X, and suppose that, for some u £ D(A), (B~ 1 Aq — X)u = r, so 



that (2.4) yields 



(2.6) -u" - fiu + u"(0)E = 7(1 -E)E. 



We will prove that (2.6), together with the fact u £ D(A), implies 7 = and u = 0, 
which shows that B~ L A — Xidx is one-to-one but not onto. Consequently, it is no 
Fredholm operator of index zero, which implies A £ <7ess(-B A , idx)- 
The general solution of (|2.6|) satisfies 



(2.7) u = c lVl +c 2V2 + U '^\ 1 E+ a2 1 i E 2 

with ci,C2 £ C, where tpi (x) = cos( v //7a;) and ip2(x) = sin(y7l:z;) if \x > 0, <px(x) = 1 
and ifi2(%) = x it (i = 0. The condition u £ L 2 (0, 00) requires c\ = c 2 — 0, and 
u(0) = tt'(Q) = then implies 

7 -u"(0) 7 7-""(0) „ 7 
a — = = 2a- 



a 2 + ix 4a 2 + fi' a 2 + fi 4a 2 + fx 

which indeed yields 7 = 0, u"(0) = 0, so that (|2.7|) provides u = 0. □ 



It is easy to generalize Theorem [2J for flow profiles which satisfy V — > c G 1R instead 
of 1. 

Corollary 2.2. Suppose that the flow profile V satisfies V —> c £ ]R, V" — » as 



2: — > 00. TTien </ie essential spectrum of (1.5) is given by 



(2.8) o- e ss(A, B) = {ll + a 2 + iaR c : ll £ [0, ex.)}. 



3. An enclosure for the set of all eigenvalues 



In this final section, we will calculate a set enclosing all eigenvalues of problem ([L5|). 
Since this set will also contain the straight line which was identified as cress (A, B) in 
Theorem 2.1, it therefore encloses the total spectrum a(A, B). 

We wish to put emphasize on the simplicity of the methods we use, which never- 
theless provide more accurate enclosures than those known in the literature Q , 0], if 
sign restrictions such as ( |l.4| ) are exploited. 

Let real constants V^, Knax, | V|max, V^' in , V^ ax be given such that 

(3.1) Knin < V(x) < V mgx , \V'(x)\ < l^'Uax, C < V" (x) < for X £ [0, co) 

for the flow profile V. Let (•, •) and || • || denote the usual inner product and norm in 
L 2 (0,oo). 
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For any eigenpair (u,X) E D(A) x C of problem (1.5) we obtain (note that B is 
positive definite) 

(Au, u) 



(3.2) 



A = 



(Bu. u) 



which is the basis of our enclosures. Partial integration and the condition u G D(A) 
H 4 (0, oo) n ff 2 (0, oo) provide 



{A 



u, u) = 



so that ( |3.2| ) yields 
(3.3) 



((-D 2 + a 2 ) 2 u, u) + iaR(V ■ (-D 2 + a 2 )u + V" ■ u, u) 
\\(-D 2 + a 2 )u\\ 2 + iaR[(Vu' ,u') - {V'u,u') + a 2 (Vu,u}] 

A = fix + iaR{p 2 - Pa) 



where 
(3.4) 

(3.5) 
(3.6) 



01 

02 



\Bu\\ 



> 



\Bu\\ 



(Bu,u) \\u\\ 
(Vu',u') +a 2 (Vu,u) 



> a" 



(u', u') + a 2 (u, u) 



(Bu,u) ' " (Bu,u) 

V(0 G [Kun, V max \ (for some C G (0, oo)), 



(V'u,u') 
(Bu, u) 



\0s\ < \V'\ 



\H\\\u'\\ 



||u'|| 2 + a 2 |M 



< 



IV'I 



2n 



From (3.3) to ( p.6[ ) we obtain the following eigenvalue enclosure result, where sums 
and products of sets are to be understood in the canonical sense. 



Theorem 3.1. All eigenvalues of problem (1.5) are contained in the set 

R 

[a 2 ,oo) + iaR[V min , V m3X ] + -\V'\ m3X ■ A, 

with A denoting the closed unit disc in C . 

Of course, separate bounds for the real and imaginary parts of the eigenvalues can 
easily be extracted (with loss of information!) from Theorem 3.1: 



Corollary 3.2. For each eigenvalue A of problem (1.5), 

a 2 -||^'|ma X < Re A 

R R 
aRV min - — | V | max < ImA < aRV max + — \V'\ max . 

Figures || and || illustrate the results of Theorem 3.1, Corollary 3.2 and Theorem 



3.3. Concretely, we have chosen V to be the Blasius profile and a = 0.179, R = 580 
here. 
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Figure 2: Illustration of Theorem 3.1 and Corollary 3.2 

The bounds given in Corollary 3.2 have been obtained in 0] already, by more involved 
analytical means. These results can be substantially improved by the following simple 
calculation, if sign restrictions such as ( [L4| ) are used: 



~({V'u,u') + {V'u,u')) 



Re (V'u, u' 

Zr 

= ~ I V"\u\*dx = 
so that fa defined in ( |3.6| ) satisfies 
(3.7) 



V'(\u\ 2 )'dx 



1 



V"{£) (u, u) (for some £ G (0, oo)) 



2 (Bu.u) 



e -i[v". V" 1 • [0 — 1 



which possibly restricts the enclosure set for /?3 used in Theorem 3.1 (which was 
obtained from ( |3.6| )). Instead of formulating a corresponding improved theorem for 
the most general situation, we concentrate on the case where V" < 0, which is true 



e.g. for the Blasius profile (see flLjp ). Then (3.7) yields Re/?3 > 0, which together 
with (|3.3|) to (|3.6|) provides the following improvement of Theorem 3.1: 



Theorem 3.3. IfV" < 0, all eigenvalues of problem ( 

R 

[a 2 ,oo) + iaR[V min , V max ] + — \V'\ 

m; 



are contained in the set 
■A-, 



where A - := {z 6 A : Imz < 0}. 

As in Corollary 3.2, we can extract separate bounds for real and imaginary parts of 
the eigenvalues, which provides here the improved bound 



aRV^ 



R 



V'\ 



; ImA < aRV m 
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In particular all eigenvalues lie "below" the essential spectrum, i.e. have imaginary 
parts less or equal aR. 



a -1/2 Rl VI 




a RV .-1/2RIVI 

mm max 



Figure 3: Illustration of Theorem 3.3 
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